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ABSTRACT 

We extend our analysis of the three-dimensional response of a vertically polytropic disk to tidal forcing 
at Lindblad resonances by including the effects of a disk atmosphere. The atmosphere is modeled 
as an isothermal layer that joins smoothly on to an underlying polytropic layer. The launched wave 
progressively enters the atmosphere as it propagates away from the resonance. The wave never propagates 
vertically, however, and the wave energy rises to a (finite) characteristic height in the atmosphere. The 
increase of wave amplitude associated with this process of wave channeling is reduced by the effect of 
the atmosphere. For waves of large azimuthal mode number m generated by giant planets embedded in 
a disk, the increase in wave amplitude is still substantial enough to be likely to dissipate the wave energy 
by shocks for even modest optical depths (t ^ 10) over a radial distance of a few times the disk thickness. 
For low-rn waves generated in circumstellar disks in binary stars, the effects of wave channeling are less 
important and the level of wave nonlinearity increases by less than a factor of 10 in going from the disk 
edge to the disk center. For circumbinary disks, the effects of wave channeling remain important, even 
for modest values of optical depth. 

Subject headings: accretion, accretion disks — hydrodynamics — waves 



1. INTRODUCTION 

Tidally generated waves can play an important role in 
the dynamics of gaseous disks and the tidal perturbers. 
The Lindblad resonances (LRs) are important sites of 
wave generation in disks (Goldreich & Tremaine 1979). 
Most previous studies have approximated the disk as two- 
dimensional (2D) by ignoring effects in the vertical direc- 
tion (perpendicular to the orbit plane). For disks whose 
structure and thermodynamic response are locally ver- 
tically isothermal, the 2D approximation accurately de- 
scribes the wave that is generated at an LR, which is in- 
deed a 2D sound wave. In many cases of astrophysical 
interest, however, the disk is thermally stratified (opti- 
cally thick) in the vertical direction. Such cases include 
the planet-forming regions of protostellar disks (Bell et al. 
1997) and cataclysmic variable (CV) disks (La Dous 1994). 

In an earlier paper (Lubow & Ogilvie 1998, hereafter 
Paper I), we analyzed the linear response of a thin, ver- 
tically polytropic disk to tidal forcing at LRs. We found 
that the dominant mode of excitation was the f mode of 
even symmetry. Close to the resonance, this mode be- 
haves compressibly and occupies the full vertical extent of 
the disk, like the 2D mode in an isothermal disk. However, 
as the wave propagates away from the resonance through 
a radial distance of order ri^/m (for resonance radius tl 
and azimuthal mode number m), it behaves incompress- 
ibly, like a surface gravity mode, and becomes confined 
close to the disk surface. The degree of confinement in- 
creases with distance from the resonance. As a conse- 
quence of this wave-channeling process, there is usually 
a strong increase in wave amplitude with radial distance 
from the LR on the scale ri^/m. Although this analysis did 
not include damping or nonlinear effects explicitly, we es- 



timated that shocks would sometimes develop and damp 
the wave. Other damping mechanisms such as radiative 
damping (Cassen & Woolum 1996) or viscous damping 
may also be of importance. 

Once the wave energy is channeled to a region close to 
the disk surface, the wave's properties are determined by 
the structure of the disk's outer layers. In particular, the 
effects of a disk atmosphere can be important. Such ef- 
fects are not well represented by a vertically polytropic 
disk. In the case of a polytrope, there is a definite surface 
on which the density and temperature drop to zero. With 
an atmosphere, the disk density drops more gradually with 
height and the temperature is approximately constant. In 
fact, there is no actual disk surface in the isothermal case. 
The density falls as a Gaussian with height above the mid- 
plane. In the limit of very high (vertical) optical depth in 
the disk, the polytropic approximation should be adequate 
because the atmosphere contains an insignificant fraction 
of the disk's mass. In the opposite limit of an optically 
thin disk, the disk structure is vertically isothermal and 
the wave launched at resonance behaves as a 2D sound 
wave. 

The purpose of this paper is to determine to conse- 
quences of a finite, non-zero optical depth on the prop- 
agation of a wave launched from an LR. For disks of op- 
tical depth of about 100, which is thought to be typical 
in planet-forming regions of protostellar disks, nearly one 
per cent of the disk's mass resides in the atmosphere. As 
a result of wave channeling, the wave energy could be con- 
tained within the atmosphere. The purpose of this paper 
is to explore the consequences of a disk atmosphere on the 
propagational properties of the waves generated at LRs. 
We aim to determine whether the wave in the disk at- 
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mosphere would propagate upwards into very low density 
regions where shock dissipation would occur. More gener- 
ally, we want to determine how much the wave amplitudes 
can be amplified by the wave-channeling process in the 
presence of a disk atmosphere as a function of disk optical 
depth. 

2. DESCRIPTION OF THE MODEL 

The semi-analytic approach taken in Paper I can be ex- 
tended to include the disk atmosphere, which is modeled as 
an isothermal layer that extends from some height above 
the mid-plane outward. Beneath this layer resides a ther- 
mally stratified polytropic layer which extends vertically 
down to the mid-plane. For convenience, we treat the tran- 
sition as abrupt, although in reality it may occur over a 
distance comparable to a scale height. 

As will be shown, physical solutions to the wave equa- 
tion in the isothermal layer can be obtained analytically. 
The upper boundary condition of the isothermal layer is 
that no waves carry energy from great distances above the 
disk mid-plane downwards toward the mid-plane. That is, 
there are no external sources of energy above the disk. The 
solutions allow for the possibility of reflection in the disk 
atmosphere and propagation to great (or infinite) heights 
in the atmosphere. 

These atmospheric solutions contain some freedom that 
needs to be constrained by the underlying polytropic layer. 
The isothermal solutions provide constraints that serve as 
an outer boundary condition for the numerical solution of 
the wave equations in the polytropic layer. Taken together, 
vertically global solutions are obtained that connect the 
two layers in a physically meaningful manner. The ratio 
of the; tcnipcratme at mid-plane to the temperature at the 
layer interface determines the disk optical depth through 
the T-T relation. 

3. EQUILIBRIUM OF THE DISK 

The disk is considered to be a steady, axisymmetric fluid 
rotating in a constant, axisymmetric gravitational poten- 
tial. Referred to cylindrical polar coordinates (r, cf), z), the 
fluid has angular velocity r2(r, z), density p(r, z), and pres- 
sure p{r,z), while the gravitational potential is ^{r,z). 
The disk is considered to be thin, so that the angular ve- 
locity can be considered independent of z and given by 



or 

while the epicyclic frequency k is given by 



^2 = 4fj2 _^ 



dr 
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The effective gravitational acceleration 
where 

9 = ^^z, 

and 0_L(r) is the vertical frequency given by 

9^2 



(1) 

(2) 

-g{r,z) e^, 
(3) 

(4) 



In particular, we will be interested in a Keplerian disk 
about a spherical mass M, for which the potential is 



$ = -GM(r2 + z'^y 



-1/2 



(5) 



and K = Q± = Q. However, the notational distinction be- 
tween il, K, and Cl±_ is retained, partly to aid the interpre- 
tation of the equations in the case of a Keplerian disk, and 
partly to allow a generalization to non-Keplerian disks. 

The local vertical equilibrium of the disk is determined 
by the equation 

d-z = 

The buoyancy frequency N{r, z) of vertical oscillations is 
given by 



N'=9 



-ypdz pdz J ' 



(7) 



where 7 is the adiabatic exponent. Throughout this pa- 
per, the self-gravitation of the fluid, its viscosity, and the 
accretion flow are neglected. 

In order to model a disk of finite optical depth, we as- 
sume that the pressure and density are related by a poly- 
tropic relation, 

p = Kp'+'/^, (8) 

in a central layer \z\ < z\ containing the mid-plane, and 
by an isothermal relation. 



p 



(9) 



in the outer layers \z\ > z\. The pressure and density are 
required to be continuous at the interface. The solution is 
then 



\z\ < zp:= po - 
\z\ > zp:^ po (^1 
for the density, and 
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for the pressure. Here po{r) and Poir) are the density and 
pressure on the mid-plane, and H = c/J7_l is the isother- 
mal scale height. The quantity Hs is the height at which 
the polytropic layer would reach zero density and pressure 
if it were not truncated, and is given by 



Hi 



2(5 + l)po 
Vl\pQ 



Furthermore, the relation 



Hi -zl = 2{s + 1)H'^ 



(14) 



(15) 



holds, because dp/dz is continuous at the interface. 

The quantity Zi/Hg is a measure of the 'weighting' as- 
signed to the polytropic layer in this mixed model. The 
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ratio of the temperature of the isothermal layer to the 
temperature at the mid-plane (for an ideal gas) is 



1 



(16) 



This may be expressed in terms of an approximate effective 
optical depth r at the mid-plane according to the relation 
(e.g. Bell et al. 1997) 



which gives 



3 ^ m 



(17) 



(18) 



The ratio Sj/Sp of the surface densities of the two lay- 
ers can be expressed in terms of zi/ Hs, but only by using 
higher transcendental functions. 

We assume that the effective adiabatic exponent 7 is 
constant and greater than unity in the central layer, but 
equal to unity in the outer layers. This reflects the fact 
that the thermal time-scale is short compared to the dy- 
namical time-scale in the isothermal layers. 

The buoyancy frequency may then be evaluated as 



\z\ < zi^.' = 



\z\>zi^.'=0, 
and is discontinuous at 



s + 1 



Z\ = Zl. 



2z' 



m - z2 



f7i 



(19) 
(20) 



4. DESCRIPTION OF THE LOCAL DISPERSION RELATION 

4.1. Basic equations 

The equations governing free, linear waves in a thin ac- 
cretion disk have been derived and solved in different cases 
by Lubow & Pringle (1993, hereafter LP), Korycansky & 
Pringle (1995), Ogilvie (1998), and in Paper L The sepa- 
ration of scales between the horizontal and vertical direc- 
tions allows all wave quantities to be expressed as WKB 
functions in r, except in the neighborhood of resonances, 
which are also turning points for the waves. In particular, 
a wave quantity X assumes the form 



X 



(r,0,z,i) ~Re{l(r,z) 



X exp 



-iujt -\- imcf) + i k{r) dr 



(21) 



where uj is the frequency eigenvalue, m is the azimuthal 
wavenumber and k{r) is the radial wavenumber. The tilde 
will be omitted hereafter. The equations satisfied by the 
Eulerian velocity perturbation {u,v,w) and pressure per- 
turbation p' at each radius may be written 



p{lu 



~,2 



k'^)u : 



dp' _ gp' 



and 



iujp' — — 7J3 iku + 



dz 

dw 
dz 



pgw, 



(22) 
(23) 

(24) 



where Vs = {'~ip/ pY^'^ is the sound speed and 

UJ = LO — mQ 



(25) 



is the intrinsic frequency of the wave. The azimuthal ve- 
locity perturbation is related to the radial one by 



where 



2Bu = 0, 



2 dr 



(26) 
(27) 



is the usual Oort parameter. 

4.2. Analytical solutions in the isothermal layer 

LP solved analytically for the vertical structure of modes 
in a purely isothermal disk by obtaining series solutions of 
the equations about z = 0. The relevant solutions were 
identified as terminating power series (i.e. polynomials) of 
either even or odd symmetry, multiplied by an exponential 
factor. In the present case, we must consider more general 
solutions because different boundary conditions apply at 
the interface between the two layers. ^ 

In the isothermal layer, equations ( |2^ ' 
rewritten as 



du 
dz 



ikw 



|) may be 
(28) 



and 



dw 
dz 



1 



(29) 



There exists a 2D mode consisting of a purely horizon- 
tal motion independent of z, which satisfies the dispersion 
relation 



For all other modes, equations (^8|) and ( |29| ) can be com- 
bined to give a single, second-order equation for w in the 
form 



dz"^ 



z 



dw 
dz 



1 



UJ 



w = Q. 



(31) 



This equation can be transformed into the parabolic cylin- 
der equation (Abramowitz & Stegun 1965, hereafter AS). 
We will be concerned with solutions for which both k and 
UJ are real. 

4.2.1. Transformation of the equation 
The substitution 



where 



z 



transforms equation ( ^T] ) into 



(32) 
(33) 

(34) 



4 



Three-dimensional waves in an accretion disk 



where 



(35) 



Two hnearly independent solutions are U{a,x) and 
V{a, x), defined by AS. The asymptotic forms as x — > +c>o 
are (AS, eq. [19.8]) 



U(a, x) - exp {~^x^) 

. 2 X 1/2 

TT 



(36) 



Via,x)^[^] x'^-i/2exp(ix2). (37) 



4.2.2. Conditions at z = ±oo 

In order to select physically acceptable solutions, we 
must examine the behavior of the wave action as \z\ — > oo. 
A suitable energy wave action can be defined, whose den- 
sity, averaged over time, is 



(38) 



and whose vertical flux, also averaged over time, is 

T'f '^ = Re {w*p') . (39) 

In the case of the isothermal layer, the limiting behavior 
of the wave action as z — s- -|-oo is 



oc x^ 

OC X 



jr(o) ^ ^-2a / 1 ^2^ 



"i-2''exp(-ix2), 
exp {-^x^) 

for the solution U{a,x), and 

^(°)oca;i+2"exp(ix2), 
oca;2'^exp(ia;2) 



(40) 
(41) 



(42) 
(43) 



for the solution V{a,x). It is clear that only the solution 
U{a,x) is acceptable at z = +oo. Similarly, only the so- 
lution U{a, —x) is acceptable at z = — oo. These solutions 
represent evanescent waves at infinity. 

4.2.3. Modes in a purely isothermal disk 

When the disk is purely isothermal without a polytropic 
layer, we recover the results of LP as follows. An accept- 
able solution for y{x) must be proportional to both U (a, x) 
and U{a,—x). However, these functions are linearly de- 



pendent only when a 
integer. In that case 



i, with n a non- negative 



U{-n - i, x) = 2-"/2i/„ (2-i/2a;) exp {-^x^) , (44) 

where Hn is the Hermite polynomial of degree n (AS, eq. 
[19.13.1]). This condition leads to the dispersion relation 



k^H^ 



(45) 



4.3. Numerical solutions in the polytropic layer 

The equations in the polytropic layer must be solved nu- 
merically, as in Korycansky & Pringle (1995). In a form 
similar to equations (Eq) and (p9), they are 



du 
dz 



1 



2z 



m 

s + l 

1 



ik 



2z2 



— K-^ _ 

w (46) 



and 



dw / 1 
dz V ik 



s + 1 \ ( (b ~ K 



m 



+ 



s + l 

7 



2z 



m-z^ 



(47) 



The boundary conditions at z = are the usual symmetry 
conditions, 

du , , 

— = w = (48) 



dz 



for an even mode, and 



dw 
dz 



(49) 



for an odd mode. At z = zi, the solutions must be 
matched on to the solutions in the isothermal layer. Both 
u and w must be continuous. 
The matching condition is 



-^ikH 






-1 r 






w 




[ )\ 





U{a,xi) 



2 



where xi — zi/H, and the prime denotes differentiation 
with respect to the (second) argument. (If this ratio di- 
verges, the matching condition is w = 0.) 

The numerical method is as follows. A value of k and a 
symmetry (either even or odd) are chosen. The values of u 
and w at z = are determined by the appropriate symme- 
try condition and by an arbitrary normalization condition. 
The value of uj must be guessed. The equations are then 
integrated from z = to z = zi , and the value of uj tuned 
so that the matching condition is satisfied. 

To evaluate the function U{a,x) and its derivative nu- 
merically, we use the following method. When a > and 
x^ + ^ 1, Darwin's expansion (AS, eq. [19.10.2]) pro- 
vides an accurate asymptotic approximation. If these con- 
ditions are not satisfied, we choose a positive integer N 
such that U{a + N, x) and U{a + N + l,x) can be eval- 
uated accurately using Darwin's expansion. These values 
are used to initialize the recurrence relation 

C/(a + n — 1, x) xU{a + n, x) + {a + n + ^)U{a + n + l,x) 

(51) 

(AS, eq. [19.6.4]), which is then iterated to determine 
U{a,x). (The recurrence is stable in this direction.) Fi- 
nally, the recurrence relation 

U'{a, x) = -^xU{a, x) - {a + ^)U{a + 1, x) (52) 



equivalent to equation (54) of LP for the case 7 = 1. 



(AS, eq. [19.6.1]) provides the derivative. 
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4.4. Asymptotic solutions in the limit kH oo 

As waves propagate radially away from the resonances 
where they are excited, the dispersion relation is often fol- 
lowed into a limit in which kH is large. This can lead to 
wave channeling, as described in Paper I, and it is impor- 
tant to determine the behavior of the modes in this limit so 
that the effects of nonlinear dissipation can be estimated. 

In a purely polytropic disk, Ogilvie (1998) showed that 
the f, p, and g modes all become trapped near the surfaces 
of the disk in this limit, and have (jj'^/ft]_ = 0{kH). In 
contrast, the r modes become trapped near the mid-plane, 
and have lo^ /Vl^ = O {{kH)~^Y (In the exceptional case 
when the disk is marginally stable to convection, the r 
modes have w^/Jl^ = O (^{kH)~^) and do not become lo- 
calized.) 

In a purely isothermal disk, the limiting behavior of the 
dispersion relation ( ^5| ) is easily found to be 

^2 ^ c2fc2 + K^ + (n+ l)nl + O {{kH)-^) (53) 

for the p modes, and 



UJ 



(n + 1) 



0{{kH) 



-4^ 



(54) 



for the r modes. Evidently the 2D mode may be considered 
to correspond to the case n = — 1 in equation (|5^). 

For the mixed model considered in this paper, it is clear 
that the limiting behavior of the r modes will agree with 
the purely polytropic case and will not be affected by the 
outer layers (except in the marginally stable case). This 
cannot be true of the other modes, however, since the sur- 
face at which they would have become concentrated is no 
longer present. Instead, the wave action moves out into the 
isothermal layer. We seek a solution in which, in accor- 
dance with the behavior of all modes other than r modes 
in the purely isothermal disk. 



c2fc2 + + (i - a)nl + O (kH) 



-2/3 



(55) 



where a is a constant to be determined. (The scaling of the 
remainder term will be seen later to be justified.) Then 



O(ikH) 



-2/3 



(56) 



The matching condition at z — zi is, from equation ( |50| ) 
ikH 



u 
w 



U'{a,xi) xi 
U{a,xi) ~ y 



O (kH) 



1/3 



(57) 



Now consider the limiting form of the equations in the 
polytropic layer. A first examination of equations (46) 
and (|3) yields 

du 

— ^ ikw 
oz 



and 



which imply 



dw 

dz 



(58) 
(59) 



dz'^ 



^2 / 



(60) 



The solution m z < zi is an evanescent WKB function, but 
there is a turning point at z = zi. In the neighborhood 
of the interface, the solution must be described using an 
Airy function. For {zi - z)/zi = O {(kH)~'^/^), one finds 



(61) 




(62) 



subject to an arbitrary normalization, and so 



w 



( 2z, y^/^ 


' Ai{Q)' 




Ai'{Q) 



= o 



at the interface. For consistency with equation (|57|), we 
require 

U'(a,xi) xi 
U{a,xi) ~ Y 



0. 



(64) 



This transcendental equation for a has infinitely many 
roots, each of which corresponds to the limiting form of 
one branch of the dispersion relation, according to equa- 
tion (H). 

5. PROPAGATION OF THE F MODE 
5.1. Numerical results 

As in Paper I, we focus of the properties of the f'' mode, 
since it is by far the dominant mode excited at an LR. In 
obtaining numerical solutions, we have considered a Ke- 
plerian disk whose polytropic layer has index s = 3 and 
adiabatic exponent 7 = 5/3. We have examined four dif- 
ferent values of the effective optical depth: r — 10, 100, 
1000 and 10000. The parameters of the four models are 
summarized in Table 1. The quantities given are the ra- 
tio Sj/Ep of the surface densities of the isothermal and 
polytropic layers, the ratio Ti/T^ of the temperature of 
the isothermal atmosphere to the temperature at the mid- 
plane, the diniensionless height zi/Hs of the interface, the 
dimensionless scale height H/Hg of the atmosphere, the 
fraction / of the torque exerted at an LR that is carried 
by the i° mode, and the eigenvalue a of equation (|^ cor- 
responding to the f° mode. 

In Figure 1 we plot the dispersion relation for the f° 
mode in the four models. This is compared with the dis- 
persion relation for a purely polytropic disk and also with 
the limiting form given by equation (|55|). It can be seen 
that the purely polytropic dispersion relation provides a 
good approximation until the wave action of the mode 
migrates into the atmosphere, at which point the approx- 
imation given by equation ( pq ) becomes good. The value 
of kHa at which this transition occurs increases with in- 
creasing r. Note that the eigenvalue a is a negative number 
whose magnitude increases with increasing r. According 
to the correspondence d = —n — i, this means that, in 
an optically thick disk, the atmospheric part of the eigen- 
function of the f'' mode does not resemble the 2D mode 
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Table 1 

Parameters of the four models. 



r 










H/H, 


/ 


a 


10 


1.266 X 10-^ 


-1 


0.7186 


0.5305 


0.2997 


0.9771 


-2.558 


100 


8.872 X 10- 


-3 


0.4041 


0.7719 


0.2248 


0.9687 


-5.751 


1000 


8.025 X 10- 


4 


0.2272 


0.8791 


0.1685 


0.9679 


-10.587 


10000 


7.670 X 10- 


5 


0.1278 


0.9339 


0.1264 


0.9678 


-18.529 



of a purely isothermal disk, but is more like the tail of a 
high-order p mode after its final node. 

We now consider the propagation of the f mode away 
from an LR, as in Paper I. Away from the resonance, the 
radial flux of angular momentum associated with the wave 
(averaged over integrated over (j) and z), 



k 



p\u\^ dz, 



(65) 



is independent of r, and this condition determines the rela- 
tive normalization of the wave at each radius. Our method 
of estimating the nonlinearity of the wave at each radius 
is to locate the value of z at which the peak of the density 
of angular momentum wave action (averaged over t), 



(66) 



occurs, and to compute the (isothermal) Mach number of 
the RMS velocity perturbation, 



M = 



HH 



1 1/2 



P/P 



(67) 



at that height. 

Since this measure depends on the overall amplitude of 
the wave, we normalize it in terms of its value at the res- 
onance as follows. In the neighborhood of the LR, the 
horizontal components of the velocity dominate, and the 
peak of the wave action density occurs on z = 0. The form 
of the inner solution is (cf. Paper I) 

u - Cu{z) [Ai{qx) ± iGi{qx)] , (68) 
V ~ iCv{z) [Ai{qx) ± iGi{qx)] , (69) 

where x — {r ~ ri^)/ri^, C is a constant, q and u{z) are 
the f-mode eigenvalue and eigenfunction discussed in Pa- 
per I, and V = —{2B/uj)u. The constant C is determined 
by asymptotic matching to the WKB solution away from 
the LR. Now the maximum value of [Ai((7a;)]^ + [Gilqx)]"^ 
is approximately 0.3603 and occurs at qx k, —1.845. This 
determines the maximum Mach number at the LR, which 
we call tMl- The scalings are such that M/Mi^ away from 
resonance is O ((iJ,/r)i/6)^ and for this reason we evaluate 
the quantity 

M \ 
.Ml/ \H, 



(70) 



As in Paper I, we have assumed that the disk is Keple- 
rian, so that the dimensionless intrinsic frequency of the 
wave is given by 



3/2 



- 1 



(71) 



where Tc is the corotation radius of the mode. We have 
also assumed that Hg ot. H oc r and E oc r~^. This vari- 
ation of parameters was selected so that the wave Mach 
number would be finite at the radial center of the disk. 
In this way, we have removed geometrical focusing effects 
from the wave-channeling effects of interest in this paper. 

In Figure 2 we plot the variation with r of the height 
Zpcak of the local vertical peak in the wave action density. 
We consider both inward propagation from the inner LR 
and outward propagation from the outer LR; we consider 
azimuthal wave numbers m — 2 and m = 10; and we 
consider the four models with different optical depths. In 
Figure 3 we plot the quantity in equation (|7^) for each of 
these cases. As Figure 2 shows, the f mode is launched at 
resonance with the peak of wave action density located at 
the mid-plane and the peak rises away from the mid-plane 
as it propagates away from the resonance. When Zpeak 
reaches the base of the atmosphere, it resides there over a 
non-zero interval in radius, as is apparent in Figure 2 by 
the flat slope of the curves. This discontinuous behavior 
of the slope is caused by our idealization of the transition 
from the polytropic interior to the isothermal atmosphere 
as occurring abruptly in the vertical direction. Notice that 
the rise of Zpeak with distance from the resonance increases 
with optical depth. However, the peak never rises to great 
heights in the atmosphere. Figure 4 shows how the wave 
action density is distributed vertically. 

We have verified that the waves are not artificially con- 
fined within the disk in our model, either by the discon- 
tinuity in the sound speed at the interface (which results 
from a change in the effective adiabatic exponent) or by 
the discontinuity in the buoyancy frequency. If the isother- 
mal layer is treated as having an adiabatic exponent equal 
to that of the polytropic layer, so that the sound speed 
is everywhere continuous, the wave migrates into the at- 
mosphere at nearly the same stage during the radial prop- 
agation. The peak of the wave action density rises only 
slightly higher than in the case of a truly isothermal at- 
mosphere. If, instead, the disk is treated as adiabatically 
stratified in both layers, so that the buoyancy frequency is 
everywhere continuous (and zero), the results are almost 
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indistinguishable from those we have presented here. We 
emphasize that radially propagating waves in an accretion 
disk are generically confined in the vertical direction be- 
cause the vertical gravity increases with height above the 
mid-plane. 

6. SUMMARY AND DISCUSSION 

We have analyzed the effects of a disk atmosphere on 
the propagation of a wave, the f mode, launched from a 
Lindblad resonance in a thermally stratified disk. The at- 
mosphere was modeled as an isothermal layer that resides 
above a polytropic region that extends down to the mid- 
plane. The polytropic region represents the (vertically) 
optically thick interior of the disk. The effects of the at- 
mosphere can be important because the f mode becomes 
progressively more confined to the disk surface with dis- 
tance from resonance. The wave energy is concentrated in 
the vertical direction as the wave propagates radially (see 
Figure 4) . The peak of wave action density in the disk fol- 
lows the trajectories shown in Figure 2 for various values 
of azimuthal wavenumber m and disk optical depth t. The 
dominant amplification of a wave occurs in the thermally 
stratified region below the disk atmosphere through the 
process of wave channeling, as was anticipated in Paper I 
(see the discussion in §8.2). The disk atmosphere lessens 
the degree of wave amplification. 

As was found in Paper I, the launched wave (f mode) 
behaves more like a surface gravity wave as it propagates 
radially away from the resonance. In a purely vertically 
polytropic disk with no atmosphere, the vertical extent of 
the surface gravity wave decreases with increasing radial 
distance from resonance. As a result, there is usually a 
substantial increase in wave amplitude on a radial scale 
of order ri^/m from resonance, through a process termed 
wave channeling. In the presence of an atmosphere, the 
effects of wave channeling are lessened (see Figure 3). The 
increase in wave amplitude becomes limited once the wave 
enters the disk atmosphere. For disks of smaller optical 
depths, the thermal stratification is small and the wave 
enters the isothermal layer on a radial scale somewhat 
smaller than ri^/m. Consequently, the wave behaves more 



like the 2D acoustic mode at most radii, in accordance 
with the expected limiting behavior of a vertically isother- 
mal disk, in which the wave is 2D at all radii and its peak 
never rises above the mid-plane. 

For disks around individual stars (circumstellar disks) 
in close binary star systems (m — 2), the increase in wave 
Mach number from resonance to disk center is typically 
less than a factor of 10, and somewhat less for lower op- 
tical depths (see Figure 3a). For disks around binary star 
systems (circumbinary disks) (m = 2) the increase in wave 
Mach number from resonance is large over a scale of order 
the resonance radius (see Figure 3c). For waves that are 
mildly nonlinear at resonance, nonlinear wave damping is 
likely to be important in such disks of modest vertical op- 
tical depth (r > 10). These effects are much stronger for 
circumbinary, rather than circumstellar disks. 

For protostellar disks (to ^ 10, see Figure 3b and 4b), 
planet-forming regions are expected to be optically thick 
(Bell et al. 1997). Typical wave velocities near a giant 
planet are likely to be mildly supersonic, with Mach num- 
ber around 0.2 (Goldreich & Tremaine 1979). In such a 
situation, wave channeling would likely lead to shocks for 
disks of moderate optical depths (of order 100) on a radial 
scale of a few times the disk thickness. 

We emphasize that the detailed nonlinear outcome of 
the wave-channeling process remains uncertain. While the 
wave action is concentrated in the polytropic layer, the mo- 
tion is approximately incompressible, but this is no longer 
true once the wave migrates into the atmosphere, where 
shocks may limit the amplitude. For this reason we have 
used the Mach number of the wave motion as a measure of 
nonlinearity. The development of shocks or other nonlin- 
ear wave phenomena under these circumstances could be 
usefully investigated using numerical simulations. 
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Fig. 1. — Local dispersion relation {solid line) for the f° mode in a Keplerian disk consisting of a polytropic layer with s = 3 and 7 = 5/3 
matched to an isothermal atmosphere. The four panels correspond to optical depths r = 10, 100, 1000, and 10000. Dotted line: dispersion 
relation for a purely polytropic disk (t = 00). Dashed line: asymptotic approximation given in equation (p5|). 

Fig. 2. — Height of the peak of the wave action density as the mode propagates away from an LR. Solid line: t = 10. Dotted line: 
T = 100. Dashed line: r = 1000. Dot-dashed line: t = 10000. Panel (a): inner LR, m = 2. Panel (b): inner LR, m = 10. Panel (c): outer 
LR, m = 2. Panel (d): outer LR, m = 10. 

Fig. 3. — RMS Mach number of the velocity perturbation at the peak of the wave action density, relative to its value at the resonance. The 
vertical scale is logarithmic. Panel (a): inner LR, m = 2. Solid line: r = 10. Dotted line: r = 100. Dashed line: t = 1000. Dot-dashed line: 
T = 10000. Panel (b): inner LR, m = 10. Panel (c): outer LR, m = 2. Panel (d): outer LR, m = 10. 

Fig. 4. — Eigenfunctions of the f° mode. The wave action density is plotted as a function of the vertical coordinate for two different values 
of the radial wavenumber. The modes are normalized so that the wave action density integrates to unity over the full vertical extent of the 
disk. Panel (a): kHs = 10. Solid line: t = 10. Dotted line: r = 100. Dashed line: r = 1000. Dot-dashed line: t = 10000. Panel (b): 
kHs = 100. 
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